Exchange interaction and correlations radically change behaviour of a quantum 

particle in a classically forbidden region 
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Exchange interaction strongly influences the long-range behaviour of localised electron orbitals 
and quantum tunneling amplitudes. It produces a power-law decay instead of the usual exponential 
decrease at large distances. For inner orbitals inside molecules decay is r~^, for macroscopic systems 
cos {k f r)r~ , where fc/ is the Fermi momentum and u = 3 for ID, i/ —3.5 for 2D and u =4 for 3D 
crystal. Correlation corrections do not change these conclusions. Slow decay increases the exchange 
interaction between localised spins and the under-barrier tunneling amplitude. The under-barrier 
transmission coefficients in solids (e.g. for point contacts) become temperature-dependent. 



PACS numbers: 31.15.xr , 71.15.-m , 71.70.Gni 



I. INTRODUCTION 

One of the first famous results of Quantum Mechan- 
ics was that a particle may tunnel through a potential 
barrier. The tunneling amplitude is exponentially small 
in the classical limit. As we will see below this result 
may be incorrect if we take into account the exchange in- 
teraction. The exchange interaction is described by the 
non-local (integration) operator, and the well-known the- 
orems proven for the Schroedinger equation with a local 
potential U{r) are violated if we add the exchange term 
(or any other non-local operator) . A similar effect is pro- 
duced by the correlation corrections. In this letter we 
consider the influence of the exchange interaction and 
correlations on an electron orbital in an atom, molecule 
or solid. The tunneling amplitude is still small in the 
classical limit, however the decay of the orbitals in the 
classically forbidden area is much slower {r^") and de- 
pends on the dimensionality of the system. 

The Hartree-Fock equation for an electron orbital ^'(r) 
in an atom, molecule or solid has the following form: 



2m dr^ 



*(r) + ([/(r) - £;)*(r) K{r) (1) 



if(r) = ^*,(r) J vl,^(r')t^^*(r')dr' 



(2) 



Here the summation runs over all electron orbitals ^'g(r) 
with the same spin projection as ^'(r). Now consider, for 
example, an inner electron atomic orbital Is. The solu- 
tion of the Schroedinger equation in potential U{r) has a 
very small range ub/Z where Z is the nuclear charge. 
Outside this range the orbital decreases exponentially 
as exp (— rZ/flf,). In the Hartree-Fock equation ^ such 
rapid decay is impossible if an atom has more than two 
electrons. Indeed, if ^'(r) ~ exp(— rZ/of,) the left-hand- 
side of eq. Il]) would be exponentially small while the 
right-hand-side is still large since K{r) in eq. ^ con- 
tains higher orbitals ^'q(r) which have larger range. The 



behaviour of the inner Hartree-Fock orbitals inside atoms 
have been studied analytically (in the semiclassical ap- 
proximation) and numerically in Ref. [l| (see also section 
IIA). The dependence on the radius r can be found from 
the multipole expansion of |l/(r — r')| in _ftr(r); the slow- 
est decay normally comes from the dipole term (~ r'/r^) 
and/or last occupied orbital 4'q(r), K ^ 5'g(r)/r^. The 
extra nodes appear since the orbitals 'l'g(r) oscillate. For 
example, the Is orbital in Cs atom has 3 nodes [1] (with- 
out the exchange term a ground state has no nodes). 
The existence of extra nodes in solutions of Hartree-Fock 
equations was also mentioned in Ref. Outside the 

atom all orbitals decay with exponential factor for an 
external electron 

Inside solids there are electrons in the conduction band 
which occupy the whole crystal. It has been pointed 
out in Ref. [ij that the exchange interaction between 
localised bound electron and the conduction band elec- 
trons leads to a power tail of the bound electron or- 
bital. The effect of the exchange interaction K{r) has 
been estimated in the free band electron approximation 
^'q(r) = exp (iq ■ r). An orbital of a bound electron de- 
creases at large distances as ill] 

^'(r) - cos(/c/r)/r^ (3) 

where kf is the Fermi momentum. Note that this solu- 
tion does not contradict to the Bloch's theorem since we 
consider localised bound electron (e.g. on an impurity 
atom) which does not belong to any electron band in the 
periodic potential. It is curious that Is orbital of an atom 
placed in a crystal has infinite number of oscillations. 

The derivation of this expression assumes the presence 
of a partly filled conduction electron band. However, in 
atoms and molecules of any length the exchange enhance- 
ment of the inner orbital tail may be mediated by a com- 
plete electron shell. The question is: can the exchange 
enhancement in solids be mediated by a nonconducting 
electron band? A special interest in this problem may be 
motivated by spintronics and solid state quantum com- 
puters based on spin qubits. The long-range tail of the 
wave function could, in principle, lead to an enhance- 
ment of the exchange spin-spin interaction between the 
distant localised spins, and enhancement of the under- 



2 



barrier tunneling amplitude. 

A special feature of the "long-tail" mechanism is that 
the state of the band electrons does not change, i.e. there 
is no need to have polarization of the conduction band by 
the localised spin. The mediating band electrons produce 
the mean exchange field K{r) in eq. ^ only. Therefore, 
this "long-tail" effect is different from other effects like 
the RKKl interaction and the double exchange spin- 
spin interaction suggested by Zener [6| (see also Refs. 
0, Q and description of Anderson and Kondo problems 
, e.g., in the book [9]). 

To investigate this problem in the present paper we 
perform calculation of the tail using the Bloch waves and 
tight-binding band electron wave functions. 



II. ATOM 

A. Exchange 

Let us first explain how the long tail appears in atoms 
The radial equation for a Hartree-Fock electronic 
orbital ^i(r) = r4>i(r) is 

^-^&2'^^U.ff-E^mr)^K,{r) (4) 



Outside the radius of the inner orbital (r > as/Z for 
Is) the energy Ei is much larger than other terms in the 
denominator of Eq. ([8|) which are of the order of i?„ 
(since the opertor in the denominator acts on ^„). In our 
example the energy of Is is \Ei\ = x 13.6 eV=4 • 10^ 
eV while the hp energy is \En\ ^ 10 eV. Therefore, we 
can approximately write 

^' Ueff-E, 2m{Ueff-E,)dr^Ueff-E, 

(9) 

The free solution in this area may be described 
by the semiclassical (WKB) approximation, 
(^f^'^'^(^r) ^ \p\~^/'^ exp{— J \p\dr/h)-, it has the usual 
range ub/Z — 0.02aB for Is. Comparison with the 
numerical solution of the Hartree-Fock equation for Is 
orbital has shown that within ~1% accuracy it is enough 
to keep the first two terms in the expansion Eq. 
beyond the classical turning point, and only one term 
at r > lOas/Z. Similar results have been obtained 
for the Dirac-Hartree-Fock orbitals which include the 
spin-orbit interaction and other single-particle rclativis- 
tic corrections [l|. Thus we see that at large distances 
^is{r) w const ^5p(r)/r2. 



B. Correlations 



Ueff^U 



(5) 



The radial exchange term can be obtained using the mul- 
tipole expansion of l/|r — r'|. Outside the radius of an 
inner orbital (e.g. in the area r > ub/Z for Is) 



E 

A;>0,r! 



Cnkb 



'nk- 



,fc+l 



(6) 



Here C„fc are the standard angular momentum dependent 
coefficients and bnk = J ^n{f')£,i{r)dr . For the multipo- 
larity k — Q the integral bnk — due to the orthogonality 
of radial wave functions with the same angular momen- 
tum. 

Now we can discuss the large distance behaviour of the 
orbital ^i(r). We will use Is orbital in Xe atom {Z — 54) 
as an example. The last occupied shells are ...5s^5p^. 
The orbital 5s does not contribute to Ki{r) since in this 
case the multipolarity of the exchange integral is fc = 
and the orthogonality condition makes bnk = 0. The 
exchange integral ls5p has fc = 1, therefore, at r ~ as 
and outside the atom Kis{r) « C5p,i65p,i^5p(r)/r^. 

The solution of Eq. ^ may be presented as [1] 



e,(r)=4/'--(r)-Her(r) 



(7) 



h^d^ 

2m dr'^ 



+ {Ueff~E,)]-'K,{r) (8) 



The effect of the correlations may be described by the 
non-local "correlation potential" S(r, r',i?) (integration 
operator) which modifies electron orbitals (see e.g. 10] ). 
The correlation potential is defined such that its aver- 
age value coincides with the correlation correction to the 
energy, 

SE, = (i\t\i) (10) 
C{r2) = ±^^ = /s(ri,r2,S,)*.(ri)dVi . (11) 



By solving the Hartree-Fock equation for the electron 
orbital including the correlation potential S, we obtain 
"Brueckner" orbitals and energies: 

l^^^w + (^w - = + (12) 

It is easy to write the correlation potential explicitly. 
In the second-order perturbation theory in the resid- 
ual interaction there are four term. The direct term 
S'^(ri,r2,Si) is given by 



V'n (''4) # (r4)V'7 (''2) V4 (r3)-0^ (ri )-0n (ra) 
r24.ri3{Ei + en ~ - e/s) 

(13) 

Note that E'' is a single-electron and energy-dependent 
operator. At large distance this term becomes the well- 
known local polarization potential ^ (see e.g. 
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so it is not interesting for us. An interesting contribution 
comes from the exchange correlation potential J^<^^'^^ 



as 



drsdr4. 



V'i (r4 ) # (r4 ) V't (i"2 (n (rs )'i/'n (rs ) 



13) 



(14) 

In this case we have the situation similar to the exchange 
interaction. Consider, for example, the correlation cor- 
rection to the Xe Is orbital (i = Is) and n = bp. The en- 
ergy of Is is large and to make an estimate we can neglect 
£5p ~ £7 ~£/3 in the denominator of the Eq. (|13p . After the 
summation over (3 and 7 we obtain the exchange correla- 
tion term ^'"'^{r) - 2e'^ / {rEis)K{v) where K{y) is the 
usual exchange term. Therefore, at large r the correla- 
tion term is suppressed in comparison with the exchange 
term by the small factor ie^ / {rEis). For Is orbital at 
r = ttB the suppression factor is The correlations 

are more important for higher orbitals where the sup- 
pression factor is 2e^ /{rEi) and \Ei \ <C \Eis\- 

We may conclude that within the perturbation the- 
ory treatment the correlations do not produce qualitative 
changes in the properties of the long-range tail. Their ef- 
fect is similar to that of the exchange, however, the decay 
is faster (extra 1/r). 



III. ID, 2D AND 3D SYSTEMS 

If we consider a molecule instead of atom, in- 
ner electron orbital will behave the same way, 
iinner{r) « coust ^valence (r)/r2. In macroscopic systems 
there is a large number of electrons occupying the valence 
band and the contribution of different valence electrons 
interfere in the exchange term in Eq.Q. This interfer- 
ence changes the long range behaviour. 

The equation for a bound electron wave function 'l'b(r) 
in a crystal contains the exchange term from Eq.Q de- 
scribing the exchange interaction of the bound electron 
with 2F mobile electrons: 



K{v)^ j g{v-v')[ 



r — r' 



-]*b(r')rfr', (15) 



5(r-r') = ^*„(r)vl/„(r') 



(16) 



Summation goes over F mobile electron states 5'n(r) with 
the same spin projection. To account for the orthogo- 
nality condition / \I'„(r')t\l/b(r')dr' = in Eq. we 

excluded the zero multipolarity term from the Coulomb 

2 22 
integrals, replacing by ^rzrp^ — In the "exact" 

2 

expression p5|) the subtracted term ^ disappears after 

the integration over r' since / \I'„(r')'''\['f,(r')cir' = 0. 

Let us start discussion of crystals from the simplest 
problem - a ID chain of N atoms separated by distance 
a. The wave function of a mobile electron can presented 



*„(r) = L-i/2e^'="-i;fe(r), 



(17) 



where Vk (r) is a periodic function in a;-direction and L = 
Na is the length of the chain. To perform the summation 
in Eq. (jl6p analytically we neglect dependence on k in 
Vk (r) . Taking the standard set of the wave vectors A:„ — 
27111/ n = 0, ±1, ibg, where F = 2g + 1, we obtain 



sin [fc/(x - a;') 
g(r - r ) = v{y)v{y ) 



(18) 



where kp = Jt^/cl and / — F/N is the band filling fac- 
tor. Now we can find the exchange term Eq The 
leading term in the multipole expansion (r' << r) of 

j^zipj — ^ ~ ^^a*" leads to the dipole approximation 
for K{r) at large distance: 



sin(fc,x)y.'cos(fc,x').(r')*.(r')rfr' 
cos(fc/a;) / x' sm{kfx')v{r')^bir')dr'll9) 



K{v) = 



It is easy to extend the problem to 2D and 3-D cases. 
In 2D case we obtain 

, , , ..JAkfR) sin (fcf i? — 7r/4) , , 
g(^r^v')^v{r)v{v')^^^ ^^3^^ (20) 

where R = r — r' and Ji is the Bessel function. In 3D 



, v(y)v(y') , ,, sinikfR). _ , 

g(r^rO^ cos {kfR)^ ' ]. (21) 



27r2i?2 



kfR 



Substituting these results into Eq (|T5)) we obtain in the 
dipole approximation that the exchange interaction term 
decays as 



K{r) 



■,{kfr)r- 



(22) 



where fc^ is the Fermi momentum and v = S for ID, 
1/ =3.5 for 2D and =4 for 3D crystal, i.e. 1/ = {5 + d)/2 
where d = 1, 2, 3 is the dimension. 

The long-range tail in Eq (|22p is due to the exchange 
interaction between bound electrons and conducting elec- 
trons which travel freely inside the crystal and may be 
found at any distance from the bound electron. As we 
have seen in section II, the perturbation theory treatment 
of the correlations does not change our conclusions. This 
is the normal metal case where the correlations are reala- 
tively weak. In this case the long-range tail of a bound 
electron orbital is the real physical phenomenon which 
should be taken into account, for example, in calculating 
tunneling amplitudes or exchange interaction between 
distant localised spins. 

Note that the expressions (I19l22p do not vanish if the 
electron band is complete. Instead they have fast oscil- 
lations if the electron Fermi momentum kf is large. This 
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conclusion looks surprising since a complete band does 
not contribute to the conductivity. One may compare 
this crystal complete band case with a molecule where 
valence electrons present on all atoms even in the ab- 
sence of the conductivity. Therefore, one may have, in 
principle, an enhanced tunneling amplitude or enhanced 
exchange interaction between distant spins (power sup- 
pression instead of exponential suppression) even in 
non-conducting materials. However, if there are strong 
electron-electron correlations (due to the strong repulsion 
between valence electrons located at the same site), they 
transform the Bloch-Hartree-Fock (conductor) state into 
the Mott insulator state where there are no free electrons 
and no long tail. 

The long-tail effect does not appear in any approach 
where the exchange interaction is replaced by an effec- 
tive potential or by a density-dependent potential. Ap- 
proximate calculations may also lead to other incorrect 
conclusions. For example, the long-range tail for a com- 
plete band case does not appear in the tight-binding ap- 
proximation for the electron wave functions. In the tight- 
binding approximation a wave function of mobile electron 



where Z > 0, / = F/N is the band filling factor and the 
last expression is obtained for I <^ N . For Z = we have 
B{F, 0) = /. Substitution of ^(r - r') from Eq. ^ 
into Eq. (fTS]) shows that if the band is partly filled, the 
tight-binding approximation leads to the same conclusion 
K{r) ~ cos {k f r)r^ . However, for the completely filled 
band / = 1 and sin (tt/Z) — 0. This means that the long- 
range exchange term vanishes in the absence of mobile 
carriers, electrons or holes. The explanation is simple: in 
the tight-binding approximation the complete band wave 
function made of the running waves Eq. (|23p is equal 
to the antisymmetrised product of the localised electron 
wave functions \I>i(r — /a). The exchange interaction with 
the localised electrons does not produce the long-range 
tail. To compare with the Bloch wave expression one may 
say that the tight-binding result for the complete band 
corresponds to K{r) ^ sin (fcyr) = for r = la. However, 
the oscillations of K{r) do not lead to vanishing of its 
effect on the wave functions - compare with the solution 
for atomic orbitals in the previous section. 



«'„(r) = N-^''^ y e*'="'''«'i(r - la), 



I 



(23) 



where 'I'i(r — la) is the one-site wave function. The sub- 
stitution of from Eq. into Eq. and summa- 
tion over n gives the following results: 

g{v - r') = XI -^(^' ^ ~ - Za)*i(r' - ma)^ (24) 



At finite temperature conducting electrons and holes 
appear. This activates the long-tail mechanism even in 
the tight-binding approximation and makes the under- 
barrier transmission coefficient temperature dependent. 
Here it may be appropriate to recall that a temperature 
dependence of the transmission coefficient has been ob- 
served near the "0.7 (2e^//i) structure" in the point con- 
tact conductance measurements . 



exp {i2TTlF/N) - 1 exp (iTrfl) 



Ar(exp(i27rZ/iV) - 1) 



ttI 



sin (tt/Z), 
(25) 
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